AD-R167  914 


UNCLASSIFIED 


PERFORMANCE  ANALVSIS  OF  A  N I CROCONPUTER-BASED 
SINOLE-LOOP  DIGITAL  CONTROL  SVSTENIU)  AIR  FORCE  HRIGHT 
AERONAUTICAL  LABS  HRIGHT-PATTERSON  AFB  OH  N  J  GAUDER 
APR  86  AFMAL-TR-85-2874  F/Q  9/5 


AD-A167  914 


AFWAL-TR- 8  5-2074 


PERFORMANCE  ANALYSIS  OF  A  MICROCOMPUTER-BASED  SINGLE-LOOP 
DIGITAL  CONTROL  SYSTEM 


M.  Gauder 

Data  Acquisition  Group 
Technology  Branch 


April  1986 


FINAL  REPORT  FOR  PERIOD  AUGUST  1982  -  AUGUST  1984 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


(§) 


DTIC 

ELECTEI 
MAY  30  1986 


D 


cr  AERO  PROPULSION  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


I 


When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said 
drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by 
implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture, 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA) 
and  is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At 
NTIS,  it  will  be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for 
publ i cation. 


Project  Engineer  Chief,  Technology  Branch 


FOR  THE  COMMANDER 


H.  T.  BUSH 

Director  '  . 

Turbine  Engine  Division 
Acu  Propuiiion  Laboratory 

"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your 

organization,  please  notify  AFWAL/POTX  _ ,  W-PAFB,  OH  45433>6563 

to  help  us  maintain  a  current  mailing  list." 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by 
security  considerations,  contractual  obligations,  or  notice  on  a  specific 
document. 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


_ At>'*l.(*Q  9'/ 


REPORT  DOCUMENTATION  PAGE 


1*.  REPORT  SECURITY  CLASSIFICATION 

Unclassified  _ 


2».  SECURITY  CLASSIFICATION  authority 

N/A  _ _ 


2t>,  OECLASSIFICATION/OOWNGRADING  SCHEDULE 

N/A  _ 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


1b.  RESTRICTIVE  MARKINGS 

N/A 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


AFWAL-TR- 85-2074 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Aero  Propulsion  Laboratory 


6c.  AOORESS  (City.  State  and  ZIP  Code) 

Air  Rjrce  Wright  Aeronautical  Laboratories  (AFSC) 
Wright-Patterson  AFB  OH  45433-6563 


b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 

dr  applicable)  Aero  Propulsion  Laboratory  (AFWAL/POT) 

AL/POTX  Air  Force  Wright  Aeronautical  Laboratories 


7b.  ADORESS  (City.  State  and  ZIP  Code) 


Wright  Patterson  AFB  OH  45433-6563 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(tf  applicable) 

AFW AL/POTX  N/A 


8c  AOORESS  ( City.  Stale  and  ZIP  Code) 

Wright  Patterson  AFB  OH  45433-6563 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 


PROJECT 

NO. 


6220 3F  3066 

ii  title  (include  Security  ckmiftca/ioniPerfonnance  Analysis  I 
of  a  Microcomputer-Based  Single-Loop  Digital 


12.  personal  authorisi  Control  System 
Michael  J.  Gauder 


13a  TYPE  OF  REPORT  13b.  time  COVEREO  14.  DATE  OF  REPORT  (Yr..  Mo..  Day) 

Final  from  _8/82 _  to  8/84  86  April 


1«.  SUPPLEMENTARY  NOTATION 

Report  contains  harmless  computer  software 


WORK  UNIT 
NO. 


15.  page  COUNT 

79 


COSATI  CODES 


FIELD 

GROUP 

09 

02 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

Microcomputer 
Analog  computer 


19.  ABSTRACT  (Continue  on  reverse  tf  necessary  and  identify  by  block  number ) 

Digital  control  system  design  has  been  popular  in  the  aerospace  and  process  control 
industries.  Advances  in  microprocessor  technologies  have  added  new  growth  to  this 
popularity.  The  use  of  microprocessor-based  systems  for  control  purposes  places  new 
demands  on  digital  control  theory.  Since  microprocessors  are  relatively  slow  digital 
machines  and  usually  have  small  work  lengths,  it  is  necessary  to  place  importance  on  the 
ffects  of  time  delays  and  amplitude  quantization  with  respect  to  the  control  system.  T 
ord  length  and  speed  of  the  signal  processing  components;  Analog-to-Digital  (A/D)  and 
Digital- to-Analoa.  CP/A)  converters  placed  additional  constraints  on  the  digital  control 
system  performance.  In  order  to  be  able  to  analyze  the- effects  of  constraints  that  aTe 
involved  in  microprocessor-based  control  system  design,  it  will  be  necessary  to  develop 
a  basic  system  and  a  set  of  procedures  that  are  modifiable  to  a  variety  of  control  con¬ 
figurations. 

(Continued  on  Reverse) 


20  oistribution/availability  of  abstract 


UNCLASSI F  IE  O/UN  LIMITE  O  5f  SAME  AS  RPT.  C  OTIC  USERS  □ 
22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Michael  J.  Gauder 

DD  FORM  1473,  83  APR  EDITION  OF  1  JAN 


21.  ABSTRACT  SECURITY  CLASSIFICATION 

I  Unclassified 


22b  TE  LEPHONE  NUMBER 
t Include  Area  Codet 

_ 513-255-3904 

EDITION  OF  1  JAN  73  IS  OBSOLETE 


22c  OFFICE  SYMBOL 

AFW AL/POTX 

SECURITY  CLASSIFICATION  OF  THIS  PAGt 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


Block  19.  (Cont) 


vThe  subject  of  this  thesis  is  the  design,  development,  and  analysis  of  a  16-bit  micro¬ 
processor  based  digital  control  system.  The  purpose  of  the  study  is  threefold: 

1.  To  show  that  a  hybrid  computer  system,  consisting  of  16-bit  single  board  microcomputer 
and  an  analog  computer,  can  be  used  effectively  for  digital  control  studies. 

2.  To  evaluate  the  frequency  response  of  the  hybrid  system. 

3.  To  identify  and  evaluate  the  error  contributors  which  can  effect  the  performance 
of  digital  control  systems. 


Acc.v.  il'-n  For 
NT’"  --,\H 


DTIC 

£CTE| 
may  3  0  1986 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


TABLE  OF  CONTENTS 


Page 


1.0  INTRODUCTION  .  1 

2.0  DESIGN  OF  THE  DIGITAL  CONTROL  SYSTEM  .  3 

2.1  Selection  of  the  Microcomputer  and  Analog  Computer.  ...  3 

2.2  Digital  Control  System  Configuration .  3 

2.3  Derivation  of  the  Control  Algorithms .  7 

2.4  Implementation  of  the  Control  Equation  on  the  Micro¬ 
computer . 7 


2.5  Simulation  of  ¥lant  Parameters  on  the  Analog  Computer  .  .  11 

2.6  Interface  Between  the  Microcomputer  and  Analog  Computer  .  13 


3.0  CONTROL  LOOP  PERFORMANCE  ANALYSIS . 16 

3.1  Response  to  a  Step  Input . 16 

3.2  Response  to  Sinusoidal  Inputs  .  19 

3.3  Frequency  Response  and  Phase  Angle  Measurements  .  19 

3.4  Improved  A/D  Conversion . 22 

3.5  Step  Responses  of  a  Tustin  Based  Controller  .  22 

3.6  Comparison  Between  the  TuBtln-Based  and  Rattan-Based 

Controllers . . . 25 

4.0  ERROR  CONTRIBUTORS  . . 28 

4.1  Analog  Computer . 28 

4.2  D/A  Biasing . 28 

f 

4.3  A/D  and  D/A  Quantization . 29 

4.4  Vord  Length . 29 

ill 


■V.VJ 


-  «  -V-V 


■  «  -  ■  -  J 


.  . - -  * 


TABLE  OF  CONTENTS  (CONTINUED) 


Page 

30 


4.5  Computation  Delay  .  .  . 

4.6  Truncation  and  Round  Off . 34 

5.0  PLANT  SIMULATION  ON  TEXAS  INSTRUMENTS  (TI)  TMS32010  DIGITAL 

SIGNAL  PROCESSING  CHIP  .  36 

'5.1  Rationale . 36 

5.2  Plant  Digitization . 37 

5.3  Control  Loop  Configuration . 37 

5.4  Implementation  of  the  Plant  on  the  Texas  Instruments  TMS 

32010  EVM  Microcomputer . 39 

5.5  Step  Response  of  the  TMS32010  Plant . 39 

5.6  State  Space  Representation  of  the  Plant  .  43 

6.0  CONCLUSIONS . 48 

APPENDICES 

A1  Program  Listing:  8-bit  A/D;  8-bit  D/A . 50 

A2  Program  Listing:  8-bit  A/D;  12-blt  D/A . 52 

A3  Program  Listing:  12-bit  A/D;  12-bit  D/A  .  54 

A4  Program  Listing:  Computational  Delay . 56 

A5  Program  Listing:  Round  Off  Routine . 58 

A6  Program  Listing:  TMS32010  Digitized  Plant  .  60 

B1  Circuit  Diagram:  Interface  Circuit  8-blt  A/D;  8-bit  D/A  .  62^ 

B2  Circuit  Diagram:  Interface  Circuit  8-blt  A/D;  12-blt  D/A.  63 

B3  Circuit  Layout  and  Diagram:  Interface  Circuit  12-bit  A/D; 

12-bit  D/A . 65 

Cl  Wiring  Diagram:  Analog  Computer  Patch  for  Uncompensated 

Plant.  .  .  . . *.  70 


iv 


TABLE  OF  CONTENTS  (CONTINUED) 


REFERENCES 


Page 

71 


LIST  OF  FIGURES 


gure 


1.  Block  diagram  of  a  sampled-data  control  loop . . 

2.  Block  diagram  of  a  digital  control  loop  with  H(s)-1  . 

3.  Block  diagram  of  a  digital  control  loop  used  for  this  study  . 

4.  Flow  chart  of  the  software  package  for  the  digital 

controller . 

5.  Analog  computer  patching  configuration . 

6.  Interface  circuitry  for  interconnecting  the  MEX68KECB  and  the 

analog  computer  . 

7.  Experimental  unit  step  response  of  the  uncompensated  plant 

with  unity  feedback . 

8.  Unit  step  responses  with  an  8-bit  A/D  converter  and  two 

different  D/A  converters  . 

9.  Experimental  sinusoidal  response  of  the  digital  control 

system  . 

10.  Frequency  response  demonstrating  improvement  of  results  with 

12-bit  D/A  converter  . 

11.  Unit  step  responses  demonstrating  Improvement  of  results  with 

12-blt  A/D  converter  . 

12.  Comparison  between  unit  step  responses  of  Rattan-  and  Tustln- 

based  controllers  with  8-bit  A/D  converter  . 

13.  Comparison  between  unit  step  responses  of  Rattan-  and  Tustln- 

based  controllers  with  12-bit  A/D  converter . 

14.  Computation  delay  results  for  Rattan-based  controller  .  .  .  . 

15.  Computation  dtflay  results  for  Tustln-based  (T-0.04  sec) 

controller  . 

16.  Computation  delay  results  for  Tustln-based  (T»0.15  sec) 

controller . < . 


17.  Unit  step  response.  Truncation  versus  round  off 


LIST  OP  FIGURES  (CONTINUED) 

Figure  *a8e 

18.  Block  diagram  of  a  MC68000/TMS32010  digital  control  loop.  .  .  38 

19.  Flow  chart  of  the  TMS32010  aoftware  package  for  the  digitized 

plant . 40 

20.  Unit  step  response  of  TMS32010  Plant  at  T-0.15  sec . 42 

21.  State  space  representation  of  the  plant  .  44 

22.  Alternative  state  space  representation  of  the  plant  .  45 


ACKNOWLEDGEMENTS 


I  wish  to  express  my  sincere  appreciation  to  Dr.  Kullp  S.  Rattan, 
Department  of  Engineering,  Wright  State  University,  for  his  contribu¬ 
tion  to  my  academic  growth  and  the  guidance  he  provided  throughout  this 
project.  I  would  also  like  to  thank  Alok  Sarwal,  a  fellow  graduate 
student,  for  his  help  at  certain  points  in  the  project.  Special  thanks 
to  Gloria  J.  Chrlsman  for  typing  the  manuscript,  Betty  J.  Baldwin  for 
proofreading,  and  to  Robert  K.  VanHook  for  producing  the  schematic 
diagrams.  Thanks  to  the  Graphics  Department  at  Wright  State  University 
for  producing  the  data  traces. 


1.0  INTRODUCTION 


Interest  In  digital  control  has  expanded  rapidly  as  a  result  of 
low  cost  16-blt  microprocessors  and  associated  support  devices  being 
Introduced.  Digital  control  Is  an  attractive  alternative  when  consider 
lng  a  control  strategy.  Therefore ,  it  is  important  that  the  capabil¬ 
ities  and  shortcomings  of  microprocessor  based  controllers  be  fully 
understood  before  they  are  put  into  service. 

It  is  well  known  to  designers  of  control  systems  that  major  dlffi- 
culties  are  found  in  mechanization  of  the  control  algorithm.  Mechani¬ 
zation  means  the  selection  of  digital  equipment,  such  as  the  Analog- 
to-Digital  (A/D)  and  Digital-to-Analog  (D/A)  converters  and  the  word 
length  of  the  computer;  the  actual  programming  of  the  algorithm;  and 
analysis  of  various  error  sources  and  the  effects  each  has  on  the  dynam¬ 
ics  of  the  controller. 

The  digital  control  configuration,  used  for  this  study,  consists 
of:  a  digital  controller.  Implemented  on  a  Motorola  MC68000  based 

microcomputer  board;  in  series  with  an  analog  plant,  simulated  on  an 
analog  computer. 

The  steps  needed  to  meet  the  stated  objectives  are: 

1.  Selection  of  the  digital  processing  components. 

2.  Generation  of  the  software  package  which  implements  (he 
control  algorithm  on  the  microcomputer  Bystem. 

3.  Simulation  of  the  analog  plant  on  an  analog  computer,  and 
interconnection  of  the  microcomputer  and  analog  computer. 


.  < 
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4.  Evaluation  of  the  performance  of  the  control  loop  for 
several  configurations. 

5.  Identification  of  items  which  cause  degradation  in  the 
performance  of  the  control  loop. 

6.  Demonstration  of  a  totally  digital  control  loop  configu¬ 
ration  where  the  plant  is  digitized  and  simulated  on  a  high  speed 
microcomputer . 

The  results  of  this  study  will  demonstrate  the  effectiveness  of 
using  a  microprocessor  based  system  for  digital  control. 
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2.0  DESIGN  OF  THE  DIGITAL  CONTROL  SYSTEM 

2.1  SELECTION  OF  THE  MICROCOMPUTER  AND  ANALOG  COMPUTER 

.The  resources  for  performing  digital  control  studies  should  be 
readily  available  and  moderate  in  cost.  The  equipment  used  for  this 
study  was  available  for  use  at  Wright  State  University.  The  Comdyna 
GP-6  and  Electronics  Associates,  Inc.,  (EAI)  TR-20  were  the  analog  com¬ 
puters  used  throughout  the  digital  control  study  effort.  The  processor 
chosen  for  the  digital  controller  was  the  Motorola  MC68000  16-bit  micro¬ 
processor.  This  chip  is  representative  of  the  many  16-blt  microproces¬ 
sors  on  the  market,  but  it  has  several  attributes  which  made  it  more 
suitable  for  the  control  study.  The  Motorola  MEX68KECB  Educational 
Circuit  Board,  a  low  cost  MC68000  based  microcomputer  board,  was  used  as 
the  digital  controller.  This  computer  board  was  purchased  for  this 
project  and  several  other  digital  control  studies  which  will  follow. 

2.2  DIGITAL  CONTROL  SYSTEM  CONFIGURATION 

A  single-loop  sampled  data  control  configuration  is  shown  in  Figure 
1.  The  primary  components  of  the  loop  are:  D^(z),  the  digital  control¬ 
ler  which  receives  and  transmits  control  data  at  sampling  instant  T; 
Gho(s),  a  zero  order  hold  device;  G(s),  the  plant  or  device  which  is  to 
be  controlled;  and  H(s),  the  feedback  element  which  takes  the  output  of 
the  plant  to  a  summing  junction  where  the  difference  between  the  set 
point  and  plant  output,  or  the  amount  of  error  remaining  to  be 
corrected,  is  fed  back  to  the  controller  input. 
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Figure  3.  Block  diagram  of  a  digital  control 


The  digital  control  loop  used  for  this  control  study  Is  similar 
to  the  configuration  of  Figure  1,  except  that  the  feedback  element, 
H(s),  was  set  equal  to  one.  Figure  2  shows  this  configuration.  Figure 
3  is  a  block  diagram  representation  of  the  control  system  as  imple¬ 
mented  . 

The  plant  G(s),  the  unity  feedback  element  H(s),  and  the  summing 
junction  were  implemented  on  an  analog  computer.  The  digital  control¬ 
ler  was  implemented  on  the  Motorola  MEX68KECB  computer  board.  The 
signal  conversion  devices,  the  A/D  and  D/A  converters,  were  part  of  an 
interface  board  which  was  developed  for  this  project. 

The  system  characteristics  are  the  following: 


G(s) 

H(s) 


6  m 

s(s^+40S+300) 

1 


2.3  DERIVATION  OF  THE  CONTROL  ALGORITHM 


(1) 

(2) 


The  pulse-transfer  function  of  the  first  order  digital  controller, 
used  with  the  control  loop,  was  obtained  using  the  computer  aided 
frequency  matching  method  of  Rattan  [1],  The  equation  for  T-0.15 
seconds  is  given  by: 


v*> 


0.154 


z-0.523 

z-0.425 


(3) 


This  control  equation  will  be  the  reference  control  algorithm  to  which 
other  algorithms  (under  evaluation)  will  be  compared. 

2.4  IMPLEMENTATION  OF  THE  CONTROL  EQUATION  ON  THE  MICROCOMPUTER 

The  digital  control  equation  D^(z)  must  be  implemented  on  a  micro¬ 
computer.  One  method  that  is  readily  adaptable  to  computer  application 
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and  the  method  chosen  for  this  study  Is  the  representation  of  Dc(z)  as  a 
difference  equation.  Equation  (3)  can  be  written  as: 

Y (z)  0. 154z-0.081  ... 

X(z)  "  z-0.425  1  ' 

Cross-multiplying  equation  (4).  multiplying  this  result  by  z  and 

solving  for  Y(z),  we  get: 

Y(z)  -  0.425z_1Y(z)+0.154X(z)-0.081z“1x(z)  (5) 

Taking  the  inverse  z-transform  of  equation  (5)  yields: 

Y(K)  -  0.425*Y(K-l)+0. 154*X(K)-0. 081*X(K-1)  (6) 


To  implement  the  first-order  difference  equation  given  in  equation 
(6)  on  a  microprocessor,  the  coefficients  have  to  be  scaled  to  a  conven¬ 
ient  base  for  ease  of  numerical  calculation.  Since  the  word  length  of 
the  MC68000  is  essentially  16  bits,  and  none  of  the  coefficients  in 
the  difference  equation  are  greater  than  one,  32767(2^)  was  chosen  as 
the  base  for  all  coefficients  to  maximize  word  length  utilization  (1 
sign  bit/13  magnitude  bits).  The  resulting  scaled  integer  coefficients 
were  then  converted  to  hexldecimal,  and  the  resulting  equation  (7)  is 
given  by : 

*  « 

Y(K)  -  3666*Y(K-1)+13B6*X(K)-0A5E*X(K-1)  (7) 

where  the  coefficients  for  equation  (7)  were  obtained  by: 


3666i6 

13B616 

0A5Ei6 


0.425*32767 

0.154*32767 


(8) 

(9) 

(10) 
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0.081*32767 


Nov  that  a  control  equation  is  In  a  form  that  can  be  Implemented 
on  the  microprocessor,  a  software  package  must  be  written  to  Instruct 
the  microprocessor  to  execute  a  sequence  of  steps  in  order  to  achieve 
the  desired  output.  The  software  package  developed  for  the  digital 
controller  consists  of  four  sections: 

1.  Initialization  section 

2.  Interrupt  servicing  and  data  input  section 

3.  Algorithm  section 

4.  Data  output  section 

The  Initialization  section  establishes  the  appropriate  configura¬ 
tion  for  the  microprocessor  and  its  support  chips.  Some  of  the  opera- 

« 

tions  performed  are:  programming  the  peripheral  interface  adapters 
(PIAs),  initializing  the  programmable  interrupt  timer  (PI/T),  and 
setting  initial  conditions  for  the  control  equation.  The  last  operation 
of  the  Initialization  section  is  to  enable  Interrupts,  enter  the  halt 
mode,  and  wait  for  an  interrupt  to  occur. 

The  interrupt  service  and  data  input  section,  which  begin  at  each 
sampling  Instant  by  acknowledging  the  Interrupt,  resets  the  interrupt 
device  and  reads  the  data  value,  X(K),  to  be  processed. 

The  algorithm  section  calculates  Y(K)  based  on  the  control  equa¬ 
tion  programmed  on  the  microprocessor,  outputs  the  results  to  the  D/A 
for  use  by  the  plant,  and  stores  appropriate  values  of  Y(K-l)  and  X(K-l) 
for  the  next  enumeration.  The  last  operation  performed  is  again  enabl¬ 
ing  interrupts,  forcing  the  processor  to  enter  the  wait  mode  until  the 
next  sampling  Instant.  The  flow  chart  of  the  software  package  is  shown 
in  Figure  4. 
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Flow  chart  of  the  software  packag 
for  the  digital  controller 


2.5  SIMULATION  OF  PLANT  PARAMETERS  ON  THE  ANALOG  COMPUTER 


The  analog  computer  provides  a  convenient  method  for  Implementing 
the  summing  junction,  the  unity  feedback  element,  and  the  plant  charac¬ 
teristics.  It  contains  a  variety  of  active  and  passive  components  which 
can  be  externally  configured  through  a  patchboard  to  simulate  the 
desired  transfer  function.  The  transfer  function  of  the  plant  can  be 
converted  Into  an  analog  computer  program  as  follows: 


r/_\  .  coo. .  6 m 

W  Y(8)  8^*408^+3008 

Cross  multiplying  equation  (11),  we  get: 

s3C(s)+40s2C(s)+300sC(b)  -  6000Y(s) 


(ID 


(12) 


Inverse  Laplace-transform  of  equation  (12)  yields: 

iiiO  +  4p  d2C<t)  +  300  d££t)  .  6W0Y(t)  (13) 

dt  dt*  dt 

In  order  to  assure  that  the  rate  of  change  of  C(t)  Is  consistent 
with  the  dynamic  properties  of  the  analog  computer  and  X-Y  plotter, 
equation  (13)  needs  to  be  "time  scaled"  before  it  can  be  Implemented  on 
the  analog  computer.  A  time  scaling  of  50  resulted  in  the  following 
equation: 


s3C(s)  -  -0.8s2C(s)-0.12sC(s)+0.048Y(s)  (14) 

This  Laplace-transform  representation  of  C(s)  can  now  be 
patched  on  the  analog  computer  using  the  configuration  shown  In  Figure 
5. 
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2.6  INTERFACE  BETWEEN  THE  MICROCOMPUTER  AND  ANALOG  COMPUTER 

Interface  circuitry,  which  would  permit  interconnection  of  the 
microcomputer  board  and  the  analog  computer,  was  developed  for  this 
study.  This  circuitry  consisted  of:  the  A/D  and  D/A  converters  and 
associated  circuitry;  two  Peripheral  Interface  Adapters  (PIA) ,  one 
programmed  as  an  input  port  (PIA1)  and  one  programmed  as  an  output  port 
(PIA2);  and  devices  used  for  chip  enable  circuits.  Figure  6  shows  a 
block  diagram  representation  of  the  interface  circuit.  Interconnection 
between  the  Interface  circuit  and  the  microcomputer  was  accomplished 
with  50  pin  ribbon  connectors  and  two  specially  made  patchcords  for 
connection  to  the  analog  computer.  The  interface  circuit  was  easily 
modifiable  for  different  A/D  and  D/A  configurations.  Let  us  take  a 
closer  look  at  each  of  the  blocks  of  Figure  6.  A  substantial  amount  of 
time  went  into  the  design  of  the  Interface  circuitry  so  a  little  more 
detailed  description  is  called  for  at  this  time. 

Motorola  MC6621  Peripheral  Interface  Adapters  (PIA)  were  used  as 

the  bus  Interface  devices  since  the  MC68000  contained  control  lines 

which  would  permit  easy  Interconnection  and  operation.  When  a  memory 

location  above  030000  Hex  was  accessed  on  the  MEX68RECB,  the  MC68000 

microprocessor  would  enter  the  synchronous  mode  of  operation.  The  Valid 

Memory  Address  (VMA*) ,  an  active-low  signal  was  used  as  on  chip  enable 

signal  for  each  PIA.  Once  a  PIA  was  selected,  a  negative-going  edge  of 

the  Enable  (E)  signal  would  cause  the.  transfer  of  data.  Three  address 

9, 

lines  (Al,  A2  and  A3)  were  used  to  select  the  proper  PIA  and  the  periph¬ 
eral  register.  Each  PIA  needed  to  be  Initialized  before  it  could  be 
used  to  transfer  data.  Writing  the  proper  data  to  the  Control  Registers 
and  Data  Direction  Registers  would  set  up  each  bit  of  the  selected  port 


as  an  input  or  output.  Bits  PA0  through  PA7  were  programmed  as  input 
lines  on  PIA1.  Bits  PB4  through  PB7  were  also  used  when  a  12-bit  A/D 
converter  was  to  be  connected  into  the  circuit.  PB0  of  PIA1  was  used  to 
provide  a  start  convert  signal  to  the  12-blt  A/D  converter.  PIA2  was 
programmed  in  a  similar  manner  as  PIA1  except  the  peripheral  ports  were 
used  as  output  lines.  The  Read/Vrlte  (R/W*)  would  determine  the  direc¬ 
tion  of  the  data  transfer.  When  the  R/W*  line  was  a  high  logic  level, 
data  was  transferred  from  the  A/D  converter,  to  a  CPU  register.  When 
this  line  was  a  low  logic  level,  data  was  transferred  from  a  CPU  regis¬ 
ter  to  the  D/A  converter. 

The  S-blt  devices  used  for  the  first  hardware  configuration  were 
National  Semiconductor  ADC0800  8-bit  successive  approximation  A/D 
converter;  the  D/A  converter  was  the  8-bit  DAC0808.  The  12-bit  D/A  con¬ 
verter  selected  for  the  second  configuration  was  the  National  Semicon¬ 
ductor  DAC1218.  The  8-bit  A/D  converter  of  the  two  previous  arrange¬ 
ments  was  replaced  with  an  Analog  Devices  AD572  12-blt  successive 
approximation  A/D  converter.  These  conversion  devices  were  selected  for 
use  since  they  were  representative  of  current  technology  and  readily 
available  for  use  in  the  laboratory. 


3.0  CONTROL  LOOP  PERFORMANCE  ANALYSIS 
3.1  RESPONSE  TO  A  STEP  INPUT 

There  were  three  signal  conversion  configurations  used  for  this 
control  study.  The  first  configuration  was  an  8-bit  A/D  converter  and 
an  8-bit  D/A  converter  arrangement.  The  second  configuration  was 
similar  to  the  first  except  that  the  8-bit  D/A  converter  was  replaced 
with  a  12-bit  D/A  converter.  The  final  arrangement  consisted  of  a 
12-bit  A/D  converter  and  a  12-bit  D/A  converter.  As  will  be  shown 
later,  the  third  arrangement  provided  the  beat  performance,  conse¬ 
quently,  a  permanent  Interface  card  was  wire-wrapped. 

The  first  step  response  to  be  measured  was  that  of  the  uncompen¬ 
sated  plant.  Figure  7  shows  the  response  obtained  when  the  uncompen¬ 
sated  plant  was  subject  to  a  unit  step  input.  Notice  that  the  step 
response  demonstrates  the  classical  overshoot  and  oscillations  associ¬ 
ated  with  an  underdamped  system.  Also  notice  that  the  steady-state 
value  of  the  plant  output  is  not  1  volt,  but  slightly  less.  This  fact 
will  be  discussed  further  in  Section  4.1. 

The  unit  step  responses  of  Figure  8  are  for  the  compensated  plant 
for  the  first-  and  second-signal  conversion  configurations.  The  over¬ 
shoot  Is  reduced  significantly,  as  would  be  expected  with  a  com¬ 
pensated  plant.  The  steady-state  oscillations  observed  will  be  discus¬ 
sed  more  thoroughly  in  Section  4.3. 


8  Bit  A/D 


3.2  RESPONSE  TO  SINUSOIDAL  INPUTS 

One  of  Che  goals  of  the  study  effort  was  to  experimentally  evaluate 
the  control  loop  response  to  sinusoidal  inputs.  This  would  provide  a 
means  of  determining  the  frequency  response  of  the  digital  control 
system.  Sinusoidal  Inputs  of  frequency  between  0  and  vere  applied 

to  the  set  point  input  of  the  control  loop.  Due  to  the  amount  of  time 
scaling  Involved  for  the  plant  simulation,  the  frequency  range  needed 
for  the  sine  waves  was  lower  than  that  obtainable  with  waveform  gener¬ 
ators  available  in  the  laboratory.  It  was  then  necessary  to  use  a 
second  analog  computer  which  generated  the  desired  sine  wave.  The 
Laplace  transform  of  the  sine  function  is  given  by: 

F(s) - y~T  (15) 

s  +u 

The  analog  computer  was  configured  for  equation  (15)  and  different 
values  for  u  were  programmed  to  yield  the  proper  input  frequency. 
Figure  9  shows  the  results  obtained  for  one  of  the  input  frequencies. 
Results  for  both  8-bit  D/A  and  22-blt  D/A  converters  are  shown  on  the 
same  plot,  along  with  the  input  frequency.  Notice  the  magnitude  atten¬ 
uation  and  phase  shift  associated  with  each  response.  The  magnitude 
attenuation  is  greater  with  the  8-bit  D/A  converter  than  with  the  12-blt 
D/A  converter.  Comparison  of  the  phase  shifts  for  each  D/A 
configuration  shows  little  difference  between  them. 

3.3  FREQUENCY  RESPONSE  AND  PHASE  ANGLE  MEASUREMENTS 

The  results  obtained  for  frequency  response  and  phase  angle  mea¬ 
surements  are  shown  in  Figure  10.  Data  for  both  D/A  converter  config¬ 
urations  are  plotted  together  with  the  theoretical  responses.  The 
theoretical  responses  for  magnitude  and  phase  angle  were  obtained  using 
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the  Interactive  control  analysis  program  TOTAL.  The  theoretical  data 
curves  provide  a  reference  to  which  the  experimental  results  can  be 
compared.  The  magnitude  plot  shows  that  the  frequency  response  with  the 
12-bit  D/A  converter  follows  the  theoretical  frequency  response  more 
closely  than  with  the  8-bit  D/A  converter.  Results  from  the  phase  angle 
plot  demonstrate  that  little  differences  exist  between  the  phase  plots 
for  the  8-bit  D/A  and  12-blt  D/A  converters,  except  at  the  highest 
frequencies  where  the  8-bit  D/A  converter  exhibited  more  deviation  from 
the  theoretical  phase  angle  curve. 

3.4  IMPROVED  A/D  CONVERSION 

All  of  the  performance  analysis  of  the  control  loop  thus  far  has 
been  with  an  8-bit  A/D  converter.  There  was  an  improvement  in  control 
system  accuracy  when  the  8-blt  D/A  converter  was  replaced  with  the 
12-bit  D/A  device.  The  8-bit  A/D  converter  will  be  replaced  with  a 
12-blt  A/D  converter,  resulting  in  the  third  control  system  configura¬ 
tion,  that  is,  12-bit  A/D  and  12-blt  D/A  converters.  The  remaining 
performance  tests  were  based  on  this  configuration.  Figure  11  contains 
step  responses,  one  where  the  loop  uses  an  8-bit  A/D  converter  and  one 
where  the  loop  uses  a  12-bit  A/D  converter.  The  step  response  associ¬ 
ated  with  the  * 12-blt  A/D  device  exhibits  slightly  less  steady-state 
oscillation  than  with  the  8-bit  A/D  converter.  The  improvement  in 
control  loop  performance  (obtained  with  Improved  A/D  conversion)  is  not 
as  pronounced  as  the  Improvement  demonstrated  with  Improved  D/A  conver¬ 
sion. 

3.3  STEP  RESPONSES  OF  A  TUSTIN  BASED  CONTROLLER 

Control  loop  performance  was  demonstrated  with  several  different 
hardware  configurations,  but  all  of  them  with  the  Rattan-based  control 
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Theoretical  Output 


algorithm.  The  Tustln  transformation  or  bilinear  transformation  [4]  as 
It  Is  commonly  known,  provides  another  means  of  obtaining  a  discrete 
system  from  the  continuous  system.  The  continuous  controller  on  which 
the  Rattan  algorithm  was  based  Is  given  by: 

CM  •  t. 322  <*» 

Substituting:  2  z  -  1 

8  T  z  ♦  1  v  n 

Into  equation  (16)  and  using  the  appropriate  value  for  the  sampling 
period  (T),  the  result  obtained  Is  a  digitized  controller  of  the  Bame 
order.  The  Tustin-based  controllers  for  T-0.15  seconds  and  T»0.04 
seconds  are  given  by  equations  (18)  and  (19),  respectively  as: 

»<*>  ■  »-352  ft :  <i8> 

■>(')  -  »-33«  a  I  Vvt)  (I9> 

The  Tustln  controller  equations  can  be  rearranged  and  the  coefficients 
converted  to  hexldecimal,  as  previously  demonstrated  with  the  Rattan 
controller,  to  obtain  the  control  algorithms: 

Y(K)  -  74B2Y(K-1)+2D04X(K)-21B7X(K-1)  (20) 

Y(K)  -  7CE3Y(K-1)+2A40X(K)-2721X(K-1)  (21) 

.  * 

for  T»0.15  seconds  and  T*0.04  seconds,  respectively.  The  Tustln  control 

algorithms  were  Implemented  on  the  digital  controller  by  changing  the 

memory  locations,  which  contained  the  associated  coefficients. 
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3.6  COMPARISON  BETWEEN  THE  TUSTIN  BASED  AND  RATTAN  BASED  CONTROLLERS 

The  step  responses  of  Figures  12  and  13  demonstrate  the  significant 
response  variations  between  the  Rattan  and  Tustin  control  algorithms. 
The  lessons  learned  from  this  indicate  that  for  a  given  hardware  config¬ 
uration,  variations  in  the  control  algorithm  can  have  a  significant 
effect  on  the  overall  performance  of  the  control  loop.  The  Tustin  con¬ 
trol  algorithms  seem  to  be  more  sensitive  to  the  size  of  the  signal  con¬ 
version  device  than  the  Rattan  algorithm.  It  is  best  to  use  the  largest 
bit  sized  conversion  device  possible  when  implementing  a  Tustin  based 
controller  to  Insure  proper  control  loop  operation. 


4.0  ERROR  CONTRIBUTORS 


4.1  ANALOG  COMPUTER 

The  transfer  function  of  equation  (1)  is  of  type  1,  which  means 
that  the  theoretical  steady-state  error  is  equal  to  zero.  However,  the 
plant,  as  implemented  on  the  analog  computer,  was  found  to  have  an  error 
of  +50  millivolts  (mV)  when  configured  with  unity  feedback  and  a  set 
point  of  1  volt,  (v),  hence,  it  was  necessary  to  establish  a  D/A  bias  at 
"digital  zero,"  which  resulted  in  a  +50  mV  output  from  the  controller. 
This  D/A  bias  would,  in  effect,  compensate  for  the  analog  computer 
error. 

4.2  D/A  BIASING 

The  use  of  a  +50  mV  D/A  bias  is  in  itself  an  induced  error,  since 
it  is  desirable  to  have  "digital  zero"  to  the  D/A  represent  a  true  value 
of  zero  volt.  Two  problems  closely  related  to  the  D/A  biasing  error  are 
over/under  D/A  biasing  and  D/A  bit  size.  If  the  bias  was  set  to  some 
value  other  than  +50  mV,  excessive  steady-state  oscillations  would 
occur.  Care  was  taken  to  insure  the  setting  of  the  proper  D/A  bias 
prior  to  any  data  collection.  Proper  setting  of  the  required  D/A 
bias  was  difficult  at  best  with  the  8-bit  D/A  converter,  but  became  less 
of  a  problem  when  the  12-bit  D/A  converter  was  used.  Establishment  of 
the  proper  D/A  bias  insured  that  the  overall  plant  response  wotld  be 


correct. 


4.3  A/D  AND  D/A  QUANTIZATION 

An  8-bit  converting  device  has  256  discrete  values,  whereas,  a  12- 
bit  converting  device  has  4096  discrete  values  associated  with  it.  For 
a  reference  voltage  range  of  10  V  (±5V),  the  resolution  for  an  8-bit  and 
a  12-bit  converters  are  39  mV  and  2.44  mV,  respectively.  Due  to  D/A 
quantization,  the  plant  output  oscillated  between  the  D/A  output  levels, 
which  drove  it  positive  or  negative.  As  the  D/A  size  was  increased,  the 
number  of  quantization  levels  also  increased,  which  resulted  in  smaller 
increments  between  the  output  levels,  therefore,  less  steady-state 
oscillation.  Similarly,  an  increase  in  A/D  bit  size  increased  the 
digital  accuracy  and  reduced  the  input  quantization  approximation  error. 

A  comparison  of  the  D/A  output  quantization  effects  can  be  seen  in 
the  unit  step  plots  of  Figure  8.  The  reduction  in  oscillation  of  the 
12-bit  configuration  is  very  evident.  The  step  response  plots  of  Figure 
11  demonstrate  that  further  improvement  in  plant  response  was  observed 
with  a  12-bit  A/D  converter,  although  this  improvement  is  not  as  signif¬ 
icant  as  seen  with  the  D/A  converter  change. 

4 . 4  WORD  LENGTH 

Another  source  of  system  error  is  the  finite  word  length  of  the 
computer.  As  seen  previously,  the  size  of  the  signal  conversion  compo¬ 
nents  has  a  significant  effect  on  the  performance  of  the  control  loop. 
The  coefficients  of  the  control  algorithm  are  scaled  values  based 
upon  a  binary  fixed-point  numerical  representation.  As  the  internal 
precision  of  the  word  length  of  the  computer  goes  up,  so  does  the 
resolution  of  the  coefficients.  This  increased  precision  propagates 
throughout  the  calculations  so  that  the  upper  word  of  the  final  computed 
value  is  a  more  accurate  representation  than  what  would  have  been 
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obtained  using  lower  precision  numerical  representation.  It  is  the 
upper  word  of  the  final  value  which  is  sent  to  the  D/A  converter.  The 
16-bit  word  length  of  the  MC68000  found  to  be  more  than  sufficient  for 
producing  acceptable  accuracy. 

4 . 5  COMPUTATION  DELAY 

The  control  algorithm  takes  a  finite  amount  of  time  to  produce  an 
output  based  upon  a  given  input.  This  delay  is  the  amount  of  time  it 
takes  to  calculate  the  control  output  at  a  given  sampling  instant  from 
an  error  input  taken  simultaneously.  The  effects  of  computation  delay 
on  control  loop  performance  may  or  may  not  be  significant.  If  the  ratio 
of  computation  delay  to  sample  rate  is  small,  then  computation  delay 
should  not  be  a  problem.  As  this  ratio  becomes  larger,  the  effects  *of 
computation  delay  on  loop  performance  should  become  apparent.  To 
experimentally  determine  computation  delay,  it  was  necessary  to  place  a 
delay  routine  in  the  control  algorithm.  The  length  of  the  delay  was 
controlled  by  a  specific  value,  placed  In  a  register,  which  was  decre¬ 
mented  until  it  was  zero.  Computation  delay  values  of  1/4T  and  1/2T 
were  used.  Figures  14  and  15  show  results  obtained  for  the  Rattan  based 
and  Tustin-based  (T-0.04  sec.)  controllers.  The  plots  for  computation 
delays  of  1/2T  show  that  during  transient  periods,  the  plant  will  tend 
to  overshoot  more  when  compared  to  the  plots  with  computation  delay  of 
1/4T.  A  comparison  of  computation  delays  of  1/16T  and  1/8T  for  the 
Tustin-based  (T-0.15  sec.)  is  given  in  Figure  16.  The  effects  of  smal¬ 
ler  computation  delays  are  more  noticable  with  this  longer  sampling 
period  than  with  the  T=0.04  sec.  controller.  This  is  an  indication  that 
the  plant  is  sensitive  to  a  fixed  amount  of  computation  delay  since 
1/16T  of  the  T«0.15  sec.  controller  is  approximately  equal  to  I/4T  of 
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the  T-0.04  sec.  controller.  As  the  plant  approaches  steady-state*  the 
effects  of  computation  delay  diminish.  For  the  control  configurations 
of  this  study,  it  appeared  that  computation  delay  did  not  effect  the 
control  loop  significantly;  however,  this  may  not  be  true  for  other 
types  of  reference  inputs. 

4.6  TRUNCATION  AND  ROUND  OFF 

Truncation  is  the  process  of  ignoring  all  bits  less  than  the  least 
significant  bit,  whereas,  round  off  is  the  process  of  selecting  a  number 
which  is  closest  to  the  unrounded  quantity.  For  example,  the  decimal 
number  1.96  will  be  truncated  to  a  value  of  1.9  and  rounded  to  2.0  for 
two  significant  digits  of  accuracy.  The  procedures  of  truncation  and 
round  off  for  binary  numbers  are  the  same.  All  of  the  control  algo¬ 
rithms  so  far  have  used  truncation  of  the  final  output  result.  The 
final  result  obtained  was  either  a  24-bit  or  a  28-bit  result,  depending 
on  the  size  of  the  A/D  converter  used.  The  most  significant  8-  or  12- 
bits  for  the  final  value  were  passed  to  the  D/A  converter,  depending  on 
the  size  of  the  D/A  device.  The  remainder  of  the  lower  significant  bits 
did  not  contribute  to  the  magnitude  of  the  final  output  value.  A 
rounding  routine  was  written  for  the  12-bit  A/D,  12-bit  D/A  control 
configuration  to  include  the  effects  on  these  lowest  bits  in  the  final 
value.  Figure  17  shows  unit  step  responses  of  the  plant;  one  with 
rounding,  one  without.  For  the  configuration  used  in  this  control 
study,  rounding  did  not  provide  significant  improvement  in  loop  perfor¬ 
mance  as  anticipated.  However,  the  step  response  of  the  control  algo¬ 
rithm  with  rounding  did  seem  to  have  a  steady-state  value  slightly 
closer  to  the  value  of  1  volt. 
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Rattan  Controller  T«  0. 1 5  sec. 


5.0  PLANT  SIMULATION  ON  THE  TEXAS  INSTRUMENTS  (TI) 

TMS 32010  DIGITAL  SIGNAL  PROCESSING  (DSP)  CHIP 

5.1  RATIONALE 

-  The  plant  for  the  control  study  thus  far  has  been  simulated  on  an 
analog  computer,  which  has  a  time  scale  factor  of  50.  This  time  scaling 
equates  to  a  sampling  period  of  7.5  seconds.  There  could  be  several 
advantages  In  replacing  the  analog  computer  with  a  digital  computer  such 
as;  more  flexibility,  elimination  of  offset  error  and  elimination  of 
time  scaling.  The  sampling  period  of  the  control  loop  whould  then  .be 
0.15  seconds  instead  of  7.5  seconds.  The  requirement  is  that  the  plant 
must  be  able  to  be  simulated  on  a  computer,  which  would  permit  proper 
operation  of  the  control  loop  at  the  desired  sampling  rate.  One 
approach  Is  to  digitize  the  plant  using  the  Tustin  transform  with  a 
sampling  period  of  0.015  seconds  and  implement  the  resulting  digital 
transfer  function  on  a  high  speed  digital  signal  processing  computer. 
The  computer  considered  for  plant  Implementation  was  the  Texas  Instru¬ 
ments  (TI)  TMS32010  Evaluation  Module  (EVM)  and  the  TI  TMS32010  Analog 
Interface  Board  (AIB).  The  EVM  bbard  is  an  evaluation  microcomputer 

board  based  upon  the  TI  TMS32010  digital  signal  processor  chip.  The  AIB 

« 

is  a  support  board  which  provides  the  necessary  12-bit  signal  conversion 

m 

so  that  thtf  EVM  board  can  be  used  for  signal  processing  applications. 
The  combination  of  these  bolrds  would  provide  everything  needed  for 
"real-time"  digital  simulation  of  the  plant. 


5.2  PLANT  DIGITIZATION 


The  transfer  function  of  the  plant  oust  be  digitized  before  it  can 
be  implemented  on  the  TMS32010  EVM.  The  Tustin  transformation!  equa¬ 
tion  (17)  must  be  substituted  into  equation  (1).  Keeping  T  unspecified 
so  that  a  general  equation  can  be  derived  and  simplifying  the  resulting 
expression  will  give  a  digitized  transfer  function  of: 


G(z) 


6000T3[z3+3z2+3z+l] 


8 [ z  3-3z 2+ 3z- 1 ]  + 1 6  OT [ z 3-z  2-z + 1 ] +600T2 [ z 3+z 2-z- 1 ) 


(22) 


If  T  is  set  equal  to  0.015  seconds,  equation  (22),  when  simplified 
becomes 


G(z) 


Y(z) 

X(z) 


0. 02025z3+0. 06075z2+0. 06075z+0. 02025 
10.535z3-26. 265z2+21.465z-5. 735 


(23) 


Equation  (23)  will  be  Implemented  in  software  on  the  TMS32010  EVM 
board. 

5.3  CONTROL  LOOP  CONFIGURATION 

The  control  loop  configuration  (using  the  TMS32010  EVM)  is  essen¬ 
tially  the  same  as  that  of  Figure  1,  except  that  now  there  is  a  digital 
plant  instead  of  analog  plant.  Figure  18  is  a  block  diagram  of  the 
control  loop  configuration  needed  for  this  portion  of  the  control  study. 
There  are  two  major  differences  between  the  control  loop  of  Figure  1  and 
the  control  loop  of  Figure  18;  the  summing  junction  is  a  difference 
amplifier  located  on  the  digital  controller  interface  card  and  the  plant 
transfer  function,  which  is  Implemented  on  the  TMS32010  EVM/A1B  combina¬ 
tion,  is  a  sampled  data  system  operating  one-tenth  of  the  controller 
sample  rate. 


5.4  IMPLEMENTATION  OF  THE  PLANT  ON  THE  TEXAS  INSTRUMENTS  TMS32010  EVM 
MICROCOMPUTER 

Equation  (23)  oust  be  converted  to  a  difference  equation  so  that 
the  plant  transfer  function  can  be  Implemented  directly  on  the  EVM. 
Solving  equation  (23)  for  Y(z)  yields  the  following: 


Y(z)  -  2.493z“IY(z)-2.037z"2Y(z)+0.544z-3Y(z)+0.00192X(z) 
♦  0.00577z_1X(z)+0.00577z_2X(z)+0.00192z“3X(z) 


(24) 


Notice  that  the  first  two  coefficients  are  larger  than  one,  which  means 
that  scaling  must  be  employed  to  obtain  functional  values  for  the 
coefficients.  The  smallest  number  that  Is  equal  to  2n  and  larger  than 
all  of  the  coefficients  is  4.  Dividing  all  coefficients  of  equation 
(24)  is  effectively  a  normalization  process.  Taking  the  inverse  Z-tran- 
sform  and  converting  the  coefficients  to  their  representative  hexl- 
decimal  values  which  results  in  the  following: 


Y(K)  -  4FCE*Y(K-l)+BEDl*Y(K-2)+1168*Y(K-3)+0010*X(K)002F*X(K-l) 

4  (25) 

+002F*X(K-2)+0010*X(K-3) 


Equation  (25)  can  now  be  programmed  directly  into  TMS32010  assembly 

language,  employing  the  same  techniques  as  used  when  the  digital  control 

equation  was  implemented  In  software.  Figure  19  is  a  flow  chart  for 

implementation  of  the  digitized  plant  transfer  function  on  the  TMS32010 
*  « 

EVM/A1B  system.  * 

5.5  STEP  RESPONSE  OF  THE  TMS32010  PLANT 

The  unit  step  response  of  the  uncompensated  digital  plant  at 
T*0.015  sec.  was  not  obtainable  for  some  unknown  reason,  so  a  search 


Into  the  possible  problems  was  conducted.  The  software  was  checked  for 


Figure  19,  Flow  chare  of  ehe  TMS  32010  package 
for  the  digitized  plant 


any  logic  or  programming  errors,  corrections  were  made,  but  the  control 
loop  still  did  not  function  properly.  Once  the  program  had  been  thor¬ 
oughly  checked,  the  next  step  was  to  try  a  slower  sampling  rate,  in  this 
case  T«0.15  sec.  The  unit  step  response  of  Figure  20  is  that  of  the 
uncompensated  digital  plant  in  closed  loop  form  with  the  slower  sampling 
rate/  The  unit  step  response  had  proven  that  the  program  was  Indeed 
working,  since  there  is  little  difference  between  this  program  and  the 
program  for  the  digitized  plant  operating  at  T-0.015  sec.  Furthermore, 
the  plant  is  undersampled  at  T-0.015  sec.,  an  indication  that  a  higher 
sampling  rate  is  definitely  needed  for  proper  plant  representation.  The 
major  difference  between  the  two  plant  programs  is  In  the  coefficients 
of  the  difference  equation.  Closer  inspection  of  equation  (24)  shows 
that  the  ratio  between  the  largest  coefficient  and  the  smallest  coef¬ 
ficient  is  approximately  1300  to  1.  This  large  of  a  coefficient  span 
was  not  represented  accurately  with  the  fixed  point  binary  numbering 
scheme.  The  use  of  coefficient  normalization  apparently  added  to  the 
problem.  By  comparison,  the  closed  loop  representation  of  the  T»0.15 
sec.  plant  required  no  coefficient  normalization  and  the  span  of  the 
coefficients  was  smaller. 

The  problem  just  discussed  becomes  worse  as  the  sampling  rate  of  a 
system  is- Increased.  Direct  implementation  of  a  difference  equation  is 
not*  feasible  particularly  when  higher  sampling  rates  are  used.  An 
alternative  method  of  implementation  that  produces  more  manageable  fixed 
k  point  coefficients  is  needed.  One  method  which  may  work  is  to  represent 
the  plant  as  a  set  of  discrete  state  equations. 
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5.6  STATE  SPACE  REPRESENTATION  OF  THE  PLANT 

The  characteristics  of  the  plant  can  be  represented  In  standard 
state  space  form  as: 


where 


x  ■  Ax  ♦  Bu 


T 

y  -  C  x 


A 


0  1  0 

0  0  1 

0  -300  -40 


B 


0 

0 

6000 


CT  -  (1  0  0] 


(26) 

(27) 


(28) 


(29) 


(30) 


The  block  diagram  of  this  system  is  shown  In  Figure  21.  An  equivalent 

T 

system  can  be  derived  by  changing  the  B  and  C  matrices  slightly.  The 
resulting  matrices  are: 


* 


B  - 


0 

0 


(31) 


CT  -  (60  0  0] 


(32) 


The  block  diagram  of  the  alternative  state  space  form  Is  shown  In  Figure 
22.  The  discrete  state  transition  equatloip  are  given  by: 


X(K  +  l)  -  ♦  X(K)  ♦  6u(K) 


(33) 


COO  -  DX(K) 


(34) 


where 

♦  -^(SI  -  A)*1  (35) 

®  "J®  *(T  ”T>  BdT  (36> 

The  alternative  state  space  representation  can  be  implemented  in  the 
following  manner.  Portions  of  Figure  22  can  be  converted  directly  into 
discrete  state  space  form.  Forming  a  system  including  only  the  first 
two  Integraters  will  give: 


'  0  1  ' 

A  -  (37) 

-300  -40. 


Solving  equation  (35)  and  equation  (36)  yields: 


(38) 


T  0.972 

0.0in 

»  - 

(39) 

L-3.346 

0.526 J 

f  0.009  ‘ 

e  - 

(40) 

L 1.115  J 

The  discrete  state  equations  of  this  system  become: 

t 

'*.(*♦  1)1  f  0.972  0.011*1  fx.Ool  P0.0091 

1  -  M  D<*>  (A1> 

*2(K  ♦  1)  L-3.346  0.526 J  U2(K>J  Ll.llSj 


•The  remaining  portion  of  the  system  to  be  implemented  is: 


Substituting  equation  (17)  into  equation  (42)  and  using  value  of  T*0.015 
sec.  will  give: 


y  -  0.45  XjOC)  ♦  0.45  X^K-l)  +  y(K-l)  (43) 

Equation  (41)  and  equation  (43)  totally  discrlbe  the  system  characteris¬ 
tics  and  can  be  implemented  in  software  without  the  coefficient  problems 
which  were  previously  discussed.  Verification  of  the  discrete  state 
space  technique  is  left  as  an  exercise  for  future  control  studies. 


6.0  CONCLUSIONS 


The  major  objective  of  this  control  study  was  to  demonstrate  the 
effectiveness  of  using  16-bit  microprocessors  for  digital  control  appli¬ 
cations.  Emphasis  was  placed  on  control  implementation  techniques  and 
error  identification  rather  than  control  algorithm  analysis.  Control 
loop  performance  was  measured  for  several  hardware  configurations  and 
several  control  algorithm  variations. 

Sources  of  error  which  effect  the  performance  of  the  control  loop 
were  identified.  Methods  were  suggested  which  would  reduce  the  error 
effects.  Quantization  error  was  the  most  troublesome  error  encountered. 
The  use  of  larger  bit-sized  converters  reduced  quantization  error  signi¬ 
ficantly.  Computation  delay  was  shown  to  introduce  a  slight  amount  of 
error  in  the  control  loop  during  transients  as  the  amount  of  delay 
increased.  Computation  delay  did  not  seem  .to  effect  the  steady-state 
behavior  of  the  control  loop.  To  Insure  proper  performance,  the 
microprocessor  must  be  able  to  execute  the  control  algorithm  well  within 
the  sampling  period  so  that  the  effects  of  the  computation  delay  will  be 
minimized. 

The  stated  objectives  of  the  control  study  were  met.  Digital  con- 

#  - 

trol  using  microprocessors  is  practical  when  considering  a  control 
strategy.  The  Increased  execution  speed  of  the  DSP  chips  will 
undoubtedly  make  these  devices  even  more  suitable  for  more  complex 
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digital  control  applications.  It  la  recommended  that  additional  control 
studies  be  performed  which  would  exploit  the  full  capabilities  of  the 
newer  DSP  chips. 
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•64,600010023 

001062 

2I9C0000 10C60 100 

MOVE.L 

•4294,600000100 

00106A 

139C00A1000 10021 

MOVE.B 

•161,600010021 

0010C2 

4E7 1 

NOP 

0010C4 

6090 

BRA. 6 

60016C2 

0010C6 

1 390606 1060 10035 

MOVE.B 

•1,600010035 

0010CE 

1390060900630005 

MOVE.B 

•15,600030605 

00  1006 

1 3900 600000 30005 

MOVE.B 

•0,600630005 

0010DE 

7E06 

MOVES. L 

•6,07 

0010E0 

5387 

6UBQ.L 

•  1,07 

0010E2 

6690 

BNE.6 

•0010E0 

0010E4 

143960630005 

MOVE.B 

•00030005, 02 

0010EA 

143906030005 

MOVE.B 

600030005,02 

001090 

0A020060 

EOR.B 

•  128,02 

00 1094 

4662 

EXT  .M 

02 

001096 

C942 

A6L.M 

•4,02 

001098 

MOVE.B 

600030001,02 

00109E 

143900030601 

MOVE.B 

600030001,02 

601164 

4266 

CLR.L 

04 

001106 

3215 

MOVE.M 

<A5> ,01 

601108 

3614 

MOVE.M 

<A4> ,03 

00 1 10A 

3662 

MOVE.M 

02 , <A4> 

00I10C 

C3D1 

MUL6.M 

<A1> ,01 

06110E 

2C0 1 

MOVE.L 

01,  D6 
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••1110 

C302 

••1112 

DCB2 

••1114 

C7D3 

•81114 

9C03 

•01118 

E384 

•0  1 1 1A 

2A84 

•01 11C 

4B1E 

001UC 

•  8C4M1B 

•01122 

4844 

•01124 

13C40M3M89 

•• 1 12A 

E884 

•0112C 

13C40003000D 

•01132 

4280 

001134 

4281 

••1134 

4282 

•01138 

4283 

•0 1 13A 

4287 

00113C 

4E73 

•0 1 13E 

066600  IB 

001142 

40E0 

MULS.H 

<A2> ,D2 

ADD.L 

02,04 

MUL6.M 

(A3) ,03 

8UB.L 

03,04 

ASL.L 

•  1,04 

MQUE.L 

04, (AS) 

BM1.6 

•00113C 

•SET 

•27,04 

6WAP.W 

04 

HOVE.B 

ASR.L 

•4,04 

MOVE.B 

O4,40003008D 

CLR.L 

00 

CLR.L 

01 

CLR.L 

02 

CLR.L 

03 

CLR.L 

07 

RTE 

BCLR 

•27,04 

BRA. 6 

•001124 
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•01000 

,4E,DI 

<coMRm°Nsb«?^‘7/2/#4> 

•01804 

4280 

CLR.L 

DO 

001004 

7204 

MOVEQ.L 

•4,01 

001008 

740 F 

MOVEQ.L 

•15, D2 

•0100A 

I3C000030803 

MOVE.B 

DO  ,400030003 

001010 

13C00003O007 

MOVE. 8 

D0,  *00030007 

001014 

13C00003000B 

MOVE .B 

D0.400030OOB 

00  10 1C 

13C00003000F 

MOVE. 8 

D0.400030O0F 

•01022 

13C00O03000 1 

MOVE . B 

00,000030001 

001028 

13C 100030003 

MOVE. 8 

Dl,  000030003 

•0102E 

13C200030005 

MOVE. 6 

02,000030005 

001034 

13C 180030007 

MOVE.B 

01,000030007 

00103A 

4400 

NOT.  8 

00 

00103C 

13C000030009 

MOVE.B 

D0  ,080030009 

001042 

13C 1088  30006 

MOVE.B 

D1.90003000B 

•01048 

13C00003000D 

MOVE.B 

00,400030000 

00104E 

13C 10003O00F 

MOVE.B 

D1.00003000F 

001054 

13FC00 80 0003O0BO 

MOVE.B 

•128,000030000 

0010SC 

13FC000000030089 

MOVE.B 

•0,800030009 

001044 

203C000E4E1C 

MOVE.L 

•937500,00 

•0I04A 

2O7C000 10025 

MOVE .L 

*45573,  A0 

001070 

0 1C80000 

MOVEP.L 

D0 ,00000 <A0) 

001074 

13FC0000000 10035 

MOVE.B 

•0,000010035 

00107C 

4280 

CLR.L 

DO 

00107E 

4281 

CLR.L 

01 

•01080 

4282 

CLR.L 

02 

001082 

4283 

CLR.L 

03 

001084 

4287 

CLR.L 

07 

•01084 

327C2000 

MOVE .W 

•8192, A1 

•0108A 

347C2002 

MOVE -W 

•8194.A2 

00108E 

347C2004 

MOVE.W 

•8194, A3 

00  1092 

387C20 10 

MOVE .W 

•8208, A4 

•01094 

3A7C2020 

MOVE.W 

•8224 ,A5 

•0109A 

32BC3444 

MOVE .W 

*13924, <Al> 

00109E 

34BC13B4 

MOVE.W 

#5044, <A2> 

0010A2 

34BC8A5E 

MOVE.W 

•2454, (A3> 

•01OA4 

4254 

CLR.W 

<A4) 

00  10A8 

4255 

CLR.W 

<A5> 

0010AA 

13FC0040000 10023 

MOVE.B 

•44,000010023 

001082 

2 1FC0000 10C40  100 

MOVE.L 

•4294,000000100 

00 10BA 

13FC00A1000 10021 

MOVE.B 

•141,000010021 

0010C2 

4E71 

NOP 

•01OC4 

40FC 

BRA. 8 

00010C2 

•010C4 

I3FC000 1000 10035 

MOVE.B 

•1,000010035 

•01OCE 

13FC0O0F000 30005 

MOVE.B 

•15,000030005 

•010D4 

13FC0 000000 30 005 

MOVE.B 

•0,008030005 

•010DE 

7C04 

MOVEQ.L 

.  *4,07 

0010E0 

53B7 

6UB0.L 

•  1 ,  D7 

•010E2 

44FC 

BNE.S 

•os  lore 

0010E4 

143900030805 

MOVE.B 

000030005,02 

•010EA 

143900030005 

MOVE.B 

000030005, D2 

00 1090 

•AO 20080 

EOR.B 

•128,02 

•0I0F4 

4882 

EXT  .W 

02 

001094 

E942 

ASL.W 

•4,02 

00 1098 

143900030001 

MOVE.B 

•00030001,02 

00 10FE 

14390003000 I 

MOVE.B 

000030001,02 

001104 

4284 

CLR.L 

04 

001104 

3215 

MOVE.W 

(A3) ,01 

001108 

3414 

MOVE.W 

<A4> ,03 

00 1 10A 

3882 

MOVE.W 

02 , <A4> 

•  •HOC 

C3D1 

MULS . w 

(Al> ,D1 

00110E 

2C0 1 

MOVE.L 

01,04 

001110 

C5D2 

MULS ,W 

<A2> ,02 

001112 

DC82 

ADD.  L 

02.04 
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001114 

C703 

001116 

?C83 

001118 

£38  6 

00111A 

2A86 

00  1 1 1C 

4B2A 

00  1 1  IE 

08C600  16 

001 122 

2E3C000 03464 

001128 

S387 

001 12A 

66FC 

00 1 12C 

4846 

00112C 

13C60003000? 

001134 

E886 

001136 

13C6000 30000 

00U3C 

4280 

00113E 

4281 

001140 

4282 

001142 

4283 

001144 

4287 

001146 

4E73 

001148 

08860018 

00114C 

60  DE 

MULS.W 

<A3) ,03 

SUB.L 

03,06 

ASL.L 

•  1,06 

MOVE.L 

06, (A5> 

BM1  .S 

•001148 

BSET 

•27,06 

MOUE.L 

•865380,07 

SUBQ.L 

•  1  ,D7 

BNE.8 

•001128 

SWAP .M 

D6 

MOUE.B 

06 ,*0003000? 

ASR.L 

•4,06 

MOUE.B 

06, *00030000 

CLR.L 

D0 

CLR.L 

01 

CLR.L 

D2 

CLR.L 

03 

CLR.L 

07 

RTE 

BCLR 

•27, D6 

BRA.S 

•00112C 
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TUTOR  J . 

1  >  MD  1000  178 | D 1 

•01000 

46FC2000 

•01064 

4280 

•01004 

7204 

•01008 

740  F 

•0100* 

13C000030003 

•01010 

I3C00603B0O7 

001014 

13C00603660B 

00101C 

13C00003000F 

•01622 

13C60603660 1 

•  0  1628 

1 3C 100030063 

00102E 

13C200030605 

•01034 

13C 100030007 

•  0t034i 

4400 

•0103C 

13C000030609 

•01042 

13C 10603000B 

•01048 

13C00063006O 

00104E 

13C 10003B00F 

•01034 

13FC00800063000D 

•010SC 

13FC000060030009 

•01064 

203C000E4E  1C 

00106* 

207C000 10025 

001070 

0 1C60000 

001074 

13FC0000000 10035 

•0 107C 

4280 

•0107E 

4281 

001080 

4282 

•01082 

4283 

•01084 

4287 

•01084 

327C2000 

•0108* 

347C2002 

•0108E 

367C2084 

•01092 

3B7C20 10 

••1696 

3*702020 

00109* 

32BC3646 

•0 109E 

34BC13B6 

•010*2 

36BC0A5E 

0010*6 

4254 

0010*8 

4253 

0010** 

13FC0040000 10023 

001082 

21FC0000 10C66 100 

00108* 

13FC00A1000 10021 

0010C2 

4E71 

0010C4 

60FC 

0010C4 

13FC00O 1000 10035 

0010CE 

13FC006F00030005 

00 1006 

13FC000600O 30005 

0010DE 

7E06 

0010E0 

3387 

0010E2 

44FC 

00I0E4 

143900030005 

0010E* 

143900630085 

0010F0 

•A6 20080 

0010F4 

4882 

0010F4 

E942 

0010FB 

143900030001 

00I0FE 

14394003008  1 

•01104 

4286 

•01106 

3215 

•01108 

3614 

•0118* 

3882 

801(0C 

C301 

•0110E 

2C01 

•01110 

C502 

•01112 

DC82 

(R0UN00FF-7/6/84  jMJG) 
MOVE.W  tie  192,  SR 


CLR.L 

D0 

MOVED. L 

04,01 

MOVEQ.L 

013,02 

MOVE.B 

D6 , •00030003 

MOVE.B 

00 ,060030007 

MOVE.B 

D8  ,*000  3000B 

MOVE.B 

06 ,000B300BF 

MOVE.B 

00,000630001 

MOVE.B 

01,000030003 

MOVE.B 

D2, 000030005 

MOVE.B 

01,000030667 

NOT  .B 

00 

MOVE.B 

00,006636609 

MOVE.B 

01 .B0603060B 

MOVE.B 

D0 ,406630600 

MOVE.B 

D 1 ,40603060F 

MOVE.B 

•128,060030000 

MOVE.B 

•0,460036069 

MOVE.L 

•937500,00 

MOVE.L 

•65573, A0 

MOVEP.L 

D0 ,00060<A0) 

MOVE.B 

•0,400010035 

CLR.L 

00 

CLR.L 

01 

CLR.L 

02 

CLR.L 

03 

CLR.L 

07 

MOVE.W 

•8192,A1 

MOVE.W 

•8194.A2 

MOVE.W 

•8196, A3 

MOVE.W 

•B268.A4 

MOVE.W 

•8224, A5 

MOVE.W 

•13926, (All 

HOUE.M 

•5046, <A2> 

MOVE.W 

•2634, < A3) 

CLR.M 

<A4> 

CLR.W 

(AS) 

MOVE .B 

•64,400010023 

MOVE , L 

•4294,006006100 

MOVE. 8 

•161,400010021 

NOP 

BRA.S 

40010C2 

MOVE .B 

•1,406010035 

MOVE .B 

•  i'5, 460036005 

move.b 

•0, 000030005 

MOVEQ.l 

.  N6, D7 

SUBQ.L 

•  1  ,D7 

BNE.6 

40010E0 

move.b 

460030005,02 

MOVE.B 

400030005,02 

EOR.B 

•  126,02  * 

EXT  .W 

02 

ASL.N 

•4,02 

MOVE.B 

00003006 1,D2 

MOVE.B 

40003606 1,D2 

CLR.L 

D6 

MOVE .W 

< A3) ,D1 

MOVE .W 

<A4> ,03 

MOVE .W 

02 ,  <A4> 

MULS.W 

(*1> ,01 

MOVE .L 

01,06 

MULS.W 

(A2) ,D2 

ADO.L 

02,06 

58 
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•01114 

C7D3 

•Bill 6 

9C83 

•01118 

E386 

•  0  1 1 1A 

2AB6 

00  1 1 1C 

662A 

00  hie 

0606000F 

•8 1 122 

67 3C 

•01124 

5255 

00 1 124 

4646 

001126 

5246 

00112A 

0BC6B006 

00112E 

13C600030009 

001134 

E666 

001134 

13C60003000O 

00113C 

4280 

00113E 

4261 

001140 

4282 

001142 

4263 

001144 

4287 

001144 

4E73 

001146 

0B06000F 

00114C 

64 1A 

00114E 

0C470000 

00 1152 

6714 

00 1154 

5355 

001150 

4646 

001158 

5344 

00115A 

08860006 

00115E 

60  CE 

001160 

68C6001B 

001144 

4646 

001166 

60  C6 

001166 

08860016 

00116C 

4646 

00114E 

408E 

MULS.W 

<A3> ,D3 

sub.l 

03,06 

ASL.L 

01,06 

MOVE .L 

06, (AS) 

BMI  ,S 

•001148 

BTST 

015,04 

BEQ.S 

•001160 

ADDQ.M 

0 1 ,  (AS) 

SMAP.M 

06 

ADDQ.W 

01,06 

BSET 

011,06 

MOVE. 6 

06, *00030009 

ASR.L 

04,06 

MOVE.B 

04, *00030000 

CLR.L 

00 

CLR.L 

D! 

CLR.L 

02 

CLR.L 

03 

CLR.L 

07 

RTE 

BTST 

015,06 

BNE.S 

•001168 

CMP.M 

00,  D7 

BEQ.S 

•001168 

SUBQ.M 

0i,(A5> 

SMAP.M 

06 

SUBQ.M 

01,06 

BCLR 

011,06 

BRA.S 

•00112E 

BSET 

027,06 

SMAP.M 

06 

BRA.S 

•B0112E 

BCLR 

027,06 

SMP.H 

06 

BRA.S 

•0B112E 
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*«  TMS320  EVM  ASSEMBLER  «» 


> 


00010 

000 

*  THIS 

PROGRAM 

1  IMPLEMENTS  THE 

UNIT  STEP  RESPONSE  OF 

00020 

000 

*  THE 

CLOSED-LOOP*  UNCOMPENSATED  PLANT.  tHE  PLANT 

00030 

000 

*  IS  DIGITIZED 

i  USING  THE  TUSTIN  TRANSFORM  WITH  T-.1S  SEC. 

00040 

000 

*  PROGRAM  WRITTEN  BY! 

MICHAEL 

J.  GAUDER  11/12/84 

00050 

000 

AORG  3 

1 

START  PROGRAM  AT  ADDR  0003 

00060 

003 

» 

00070 

003 

*  SYSTEM  EQUATES  AND 

DATA 

ASSIGNMENTS 

00080 

003 

* 

00090 

003 

0000 

ZERO 

DATA 

>0000 

ZERO 

00100 

004 

0001 

ONE 

DATA 

>0001 

ONE 

00110 

005 

0004 

FOUR 

DATA 

>0004 

FOUR 

00120 

006 

0008 

MODE 

DATA 

>0008 

SAMPLE  DELAY ICONT.  CONVERSION 

00130 

007 

0005 

SAMRT 

DATA 

>0005 

SAMPLE  PERIOD  0.010  SEC  APPR. 

00140 

008 

7FF0 

MASK1 

DATA 

>7FF0 

MASK  FOR  INPUT  DATA 

00150 

009 

8000 

MASK2 

DATA 

>8000 

MASK  FOR  OUTPUT  DATA 

00160 

OOA 

PAO 

EQU 

0 

PORT  ADDR  FOR  AIB  CONTROL  PORT 

00170 

OOA 

PA1 

EQU 

1 

PORT  ADDR  FOR  SAMPLE  RATE  PORT 

00180 

OOA 

PA2 

EQU 

2 

PORT  ADDR  FOR  A/D  AND  D/A 

00190 

OOA 

* 

00200 

OOA 

*  COEFFICIENTS  FOR  DIFFERENCE  EQUATION 

00210 

OOA 

* 

00220 

OOA 

CCEE 

AO 

DATA 

>CCEE 

FP  VALUE  FOR  NORMALIZED  AO 

00230 

008 

A3F9 

A1 

DATA 

>A3F9 

FP  VALUE  FOR  NORMALIZED  A1 

00240 

OOC 

D387 

A2 

DATA 

>D3B7 

FP  VALUE  FOR  NORMALIZED  A2 

00250 

OOD 

2746 

A3A6 

DATA 

>2746 

FP  VALUE  FOR  NORMALIZED  A3SA6 

00260 

OOE 

7645 

A4AS 

DATA 

>7645 

FP  VALUE  FOR  NORMALIZED  A4SA5 

00270 

OOF 

« 

00280 

OOF 

*  VARIABLE  STORAGE 

00290 

OOF 

« 

00300 

OOF 

**  Y 

VALUES 

** 

00310 

OOF 

TEMP 

BSS 

1 

TEMP  YtK)  LOC. 

00320 

010 

Y1 

BSS 

1 

Y(K-l) 

00330 

Oil 

Y2 

BSS 

1 

Y(K-2> 

00340 

012 

Y3 

BSS 

1 

YCK-3> 

00350 

013 

••  X 

VALUES 

»* 

00360 

013 

XO 

BSS 

1 

X<K> 

00370 

014 

XI 

BSS 

1 

X<K-1> 

00380 

015 

X2 

BSS 

1 

X(K-2) 

00390 

016 

X3 

BSS 

1 

XCK-3> 

00400 

017 

* 

00410 

017 

*  INITIALIZE 

THE  VARIABLE 

STORAGE  TO  ZERO 

00420 

017 

* 

00430 

017 

6E00 

LDPK 

0 

LOAD  DATA  POINTER  WITH  0 

00440 

018 

7007 

LARK 

0*7 

LOAD  AUX.  REG  0  WITH  COUNT 

00450 

019 

7E03 

LACK 

ZERO 

LOAD  STARTING  VARIABLE  ADDR 

00460 

01A 

6703 

TBLR 

ZERO 

INIT  D.M. 

00470 

018 

710F 

LARK 

I  *  TEMP 

START  OF  DATA  STORAGE 

00480 

01C 

6881 

CVS 

LARP 

1 

AUX.  REG.  POINTER  SET  TO  1 

00490 

©ID 

67A0 

TBLR 

t+*0 

CLEAR  MEM.  LOC. ! INC.  ADDR»ARP»l 

00500 

01E 

F400 

BANZ 

CVS 

CONT.  UNTIL  ARO  -  0 

OIF 

001 C 

• 

00510 

020 

• 

00520 

020 

»  INITIALIZE 

DATA  MEMORY 

WITH  CONTROL  VALUES 

00530 

020 

* 

00540 

020 

7E04 

LACK 

ONE 

LOAD  ACC  WITH  ADDR  OF  ONE 

00550 

021 

6704 

TBLR 

OrtE 

STORE  MM  INTO  DM 

00560 

022 

7009 

LARK 

0*9 

LOAD  AUX.  REG  0  WITH  COUNT 

00570 

023 

7105 

LARK 

1 *FOUR 

AUX.  HAS  STARTING  D.M.  ADDR 

00580 

024 

7E05 

LACK 

FOUR 

ACC  HAS  START  ADDR  FOR  INIT 

00590 

025 

6881 

CDI 

LARP 

1 

AUX.  REG.  POINTER  -  1 

00600 

026 

67A0 

TBLR 

•4*0 

HOVE  PM  TO  DM! INC  DM  POINTER 
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00610 

027 

0004 

AOO 

ONE 

INC  PM  POINTER 

00620 

028 

F400 

BANZ 

CDI 

CONT.  UNTIL  DONE 

029 

0025 

00630 

02A 

* 

00640 

02A 

» 

ANALOG  80AR0  INITIALIZATION 

00650 

02A 

« 

00660 

02A 

700F 

LARK 

Of  >000F 

LOAD  ARO  WITH  COUNT  FOR  81 

00670 

02B 

4907 

OUT 

SAMRTfPAl 

SET  UP  SAMPLE  CLOCK 

00600 

02C 

4806 

OUT 

MODEfPAO 

8ET  UP  MODE!  START  CLOCK 

00690 

020 

* 

00700 

020 

» 

WAIT 

FOR  OATA  TO  BE  CONVERTED 

00710 

020 

» 

00720 

020 

F600 

WAIT 

BIOZ 

CNTDWN 

START  COUNTDOWN  WHEN  BIO  • 

02E 

FFFF 

00730 

02F 

F900 

B 

WAIT 

WAIT  UNTIL  BIO  -  0  BEFORE 

030 

0020 

00740 

031 

6880 

CNTDWN 

LARP 

0 

LOAD  AR  POINTER  WITH  ZERO 

007S0 

032 

F400 

BANZ 

WAIT 

IF  ARO  NE  ZERO f GO  BACK  TO 

033 

0020 

00760 

034 

700F 

LARK 

0t>000F 

RESET  COUNTDOWH  REGISTER 

00770 

035 

t 

00780 

035 

« 

REAO 

OATA 

AND  CONVERT 

TO  TRUE  BINARY 

00790 

035 

« 

00800 

035 

7F89 

CALC 

ZAC 

CLEAR  ACCUMULATOR 

00810 

036 

4213 

IN 

X0>PA2 

INPUT  DATA  FROM  A/D 

00820 

037 

2013 

LAC 

XO 

LOAD  ACC  FROM  X(K> 

00830 

038 

7808 

XOR 

HASK1 

COMPLEMENT  SIGN  BIT 

00840 

039 

5013 

SACL 

XO 

STORE  TRUE  BIN  AT  X(K> 

00850 

03A 

7F89 

ZAC 

CLEAR  ACC. 

00860 

03P 

« 

00870 

033 

« 

CALCULATE 

DIFFERENCE 

EQUATION 

00880 

03B 

* 

00890 

03B 

6A12 

LT 

Y3 
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NOTES 

X  -  DENOTES  VIKE  WRAP  PINS  FOR  BYPASS  CAPACITORS 
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50  PIN  RIBBON  CONNECTOR 
ON  WIRE  WRAP  INTERFACE  BOARD 
TO  MEX68KECB  CPU  BOARD 
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